Development of the 12-Base Short Dimeric Myogenetic Oligodeoxynucleotide That Induces Myogenic Differentiation

A myogenetic oligodeoxynucleotide (myoDN), iSN04 (5′-AGA TTA GGG TGA GGG TGA-3′), is a single-stranded 18-base telomeric DNA that serves as an anti-nucleolin aptamer and induces myogenic differentiation, which is expected to be a nucleic acid drug for the prevention of disease-associated muscle wasting. To improve the drug efficacy and synthesis cost of myoDN, shortening the sequence while maintaining its structure-based function is a major challenge. Here, we report the novel 12-base non-telomeric myoDN, iMyo01 (5′-TTG GGT GGG GAA-3′), which has comparable myogenic activity to iSN04. iMyo01 as well as iSN04 promoted myotube formation of primary-cultured human myoblasts with upregulation of myogenic gene expression. Both iMyo01 and iSN04 interacted with nucleolin, but iMyo01 did not bind to berberine, the isoquinoline alkaloid that stabilizes iSN04. Nuclear magnetic resonance revealed that iMyo01 forms a G-quadruplex structure despite its short sequence. Native polyacrylamide gel electrophoresis and a computational molecular dynamics simulation indicated that iMyo01 forms a homodimer to generate a G-quadruplex. These results provide new insights into the aptamer truncation technology that preserves aptamer conformation and bioactivity for the development of efficient nucleic acid drugs.


Introduction
Myogenic progenitor cells, called myoblasts, differentiate into myocytes, fuse into multinucleated myotubes, and eventually form myofibers during skeletal muscle development and regeneration.Myogenic differentiation is thus an initial key step for muscle growth and maintenance [1].However, aging, chronic kidney disease, cancer, or diabetes mellitus decrease the differentiation capacity of myoblasts [2][3][4][5], which is one of the reasons BioTech 2024, 13, 11. https://doi.org/10.3390/biotech13020011https://www.mdpi.com/journal/biotechfor age-and disease-associated muscle loss such as sarcopenia and cancer cachexia.Therefore, molecules that enhance myogenesis may be promising drugs to prevent the muscle wasting that is increasing in aging societies.Oligodeoxynucleotides (ODNs) have been studied as nucleic acid aptamers that bind to their targets in a conformation-dependent manner with high specificity and affinity, similar to antibodies.Nowadays, aptamers are one of the potential next-generation drugs [6].We have recently reported a series of myogenetic ODNs (myoDNs), which are single-stranded 18-base telomeric phosphorothioated (PS)-ODNs, that robustly induce myogenic differentiation of myoblasts [7] and rhabdomyosarcoma [8].
The nucleic acids containing telomeric hexamers (5 ′ -TTA GGG-3 ′ ), such as myoDNs, are known to form guanine (G)-quartet/quadruplex structures and to recruit DNA/RNAbinding proteins [6].Indeed, one of the myoDNs, iSN04 (5 ′ -AGA TTA GGG TGA GGG TGA-3 ′ ), serves as an aptamer that physically interacts with nucleolin, a multifunctional RNA-binding phosphoprotein localized in myoblast nuclei [7,9].Antagonizing nucleolin with iSN04 reverses nucleolin-inhibited translation of p53 mRNA and improves p53 protein levels, resulting in enhanced myogenesis [7].Since iSN04 restores myoblast differentiation deteriorated by diabetes mellitus [10] or cancer-secreted factors [11], it is expected to be a drug seed to reverse muscle loss in these diseases.In addition, nucleolin inhibition by iSN04 suppresses inflammatory responses [12] and induces myocardial differentiation of pluripotent stem cells [13] by modulating the β-catenin signaling pathway.These results suggest that iSN04 may be applicable to inflammatory muscle diseases such as dystrophy and to cardiac regeneration therapy.iSN04 is incorporated into myoblasts without a carrier [7].In general, single-stranded PS-ODNs can be spontaneously incorporated into the cytoplasm by an endocytic process termed gymnosis, in part due to their lower molecular weights compared to doublestranded nucleotides [14].Since shorter PS-ODNs are more favorable for intracellular uptake, sequence shortening is important pharmaceutically to improve drug efficacy and is also beneficial industrially to reduce synthesis costs [15].In this study, we attempted to shorten iSN04 while maintaining its myogenetic activity.In the previous report, computational simulations and mutational experiments revealed that the guanine stack (G-stack; GGG at the 13-15th bases) within iSN04 is the core structure.Interestingly, an isoquinoline alkaloid, berberine, binds to the G-stack and enhances the activity of iSN04 [7,9].By forming a complex, berberine is thought to shift the conformation of iSN04 to a more stable and active form, suggesting that the ODNs that configure the G-stack may serve as myoDNs.Here, we investigated the myogenetic activities of the 11~16-base non-telomeric PS-ODNs predicted to form the G-stack structure.Several PS-ODNs successfully induced human, mouse, and chicken myoblast differentiation as well as iSN04 does.These results provide new insights into the design of shorter aptamers.

Oligonucleotides
PS-ODNs in which all phosphodiester bonds were phosphorothioated to enhance nuclease resistance in cell culture were synthesized and HPLC-purified (GeneDesign, Osaka, Japan) [7,16].Unmodified oligonucleotides (DNA and RNA) without phosphorothioate for nuclear magnetic resonance (NMR) and native polyacrylamide gel electrophoresis (PAGE) were synthesized and HPLC-purified (Hokkaido System Science, Sapporo, Japan).The oligonucleotide sequences are listed in Table 1.iMyo01-iMyo08 are the myoDN candidates designed and evaluated in this study.They have two G groups (GGG or GGGG) separated by a mini linker (T or TGA), thymines at the 5 ′ end (TT or TTT), and adenines at the 3 ′ end (AA or AAA).iSN04 is a positive control that induces myogenesis.iSN14 and iSN45 are negative controls that have no myogenic capacity [7].iMyo01-CS is a reverse complementary strand DNA to iMyo01.In this study, 12nt-Ctrl is a negative control DNA that does not form a G-quadruplex, and 12nt-Ctrl-CS is a reverse complementary strand RNA to 12nt-Ctrl.All oligonucleotides were dissolved in endotoxin-free water, and an equal volume of water without oligonucleotide was used as a control.

Myoblast Culture
All myoblasts were cultured on the dishes/plates coated with collagen type I-C (Cellmatrix; Nitta Gelatin, Osaka, Japan) at 37 • C under 5% CO 2 throughout the experiments.

Immunocytochemistry
After 48 h of PS-ODN treatment, myoblasts were subjected to immunostaining for myosin heavy chain (MHC).The myoblasts were fixed with 2% paraformaldehyde, permeabilized with 0.2% Triton X-100, and immunostained with 0.5 µg/mL mouse monoclonal anti-MHC antibody (MF20; R&D Systems, Minneapolis, MN, USA).As a secondary antibody, 0.1 µg/mL of Alexa Fluor 488-conjugated donkey polyclonal anti-mouse IgG antibody (Jackson ImmunoResearch, West Grove, PA, USA) was used.Cell nuclei were stained with DAPI (Nacalai).The screening system for analyzing MHC images has been described previously [7].Briefly, fluorescence images were acquired automatically using CellInsight NXT (Thermo Fisher Scientific) and quantified using HCS Studio: Cellomics Scan software version 6.4.1 (Thermo Fisher Scientific) as follows: For hMBs and chMBs, MHC signal intensity was defined as the total MHC signal intensity divided by the total number of nuclei.For mMBs, the ratio of MHC + cells was defined as the number of nuclei in the MHC + cells divided by the number of nuclei.High-resolution MHC images (1280 × 960 pixels, 8 bit) were obtained using EVOS FL Auto microscope (AMAFD1000; Thermo Fisher Scientific) with ×10 objective magnification.Exposure time and bandpass filter (excitation/emission) were 100 ms and 470/510 nm for MHC and 15 ms and 357/447 nm for DAPI, respectively.The ratio of MHC + cells was defined as the number of nuclei in all MHC + cells divided by the total number of nuclei, and the fusion index was defined as the number of nuclei in multinucleated MHC + myotubes divided by the total number of nuclei, which were calculated using ImageJ software version 1.52a (National Institutes of Health, Bethesda, MD, USA).

qPCR
After 24 h of PS-ODN treatment, the total RNA was isolated from hMBs using Nucleospin RNA Plus (Macherey-Nagel, Düren, Germany) and reverse transcribed using a ReverTra Ace qPCR RT Master Mix (TOYOBO, Osaka, Japan).qPCR was performed using a GoTaq qPCR Master Mix (Promega, Madison, WI, USA) with a StepOne Real-Time PCR System (Thermo Fisher Scientific).The amount of each transcript was normalized to that of the tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta gene (YHHAZ).The primer sequences were described previously [7,11].Results are expressed as fold-change.

Agarose Gel Electrophoresis
An amount of 0.8 nmol PS-ODNs and 0.8 nmol berberine hydrochloride (Nacalai) were mixed in 16 µL Ham's F10 containing 144.1 mM Na + , 0.6 mM Mg 2+ , 5 mM K + , 0.3 mM Ca 2+ , 1 µM Fe 2+ , 1 µM Cu 2+ , and 0.1 µM Zn 2+ cations.The mixtures were placed at 4 • C overnight and then subjected to electrophoresis using a TAE-buffered 3% agarose gel with 0.5 µg/mL ethidium bromide (EtBr).For monochromatic images, the gels were illuminated with 365 nm ultraviolet (UV) light and the images were captured using ImageQuant LAS 500 with a 560 nm emission bandpass filter (GE Healthcare).For colored images, the gels were illuminated by 302 nm UV light and the images were taken by a digital still camera without any filters to detect the 530 nm yellow emission from berberine and the 620 nm red emission from EtBr [7].

NMR
iMyo01-DNA was annealed by heating at 95 • C for 5 min followed by snap-cooling on ice.Annealed iMyo01 (final concentration of 0.7 mM) was dissolved in 20 mM sodium phosphate (pH 6.5) using a centrifuge concentrator (Vivaspin 3000 MWCO; Sartorius, Göttingen, Germany).For the K + condition, concentrated KCl was added to the sample to a final concentration of 50 mM.NMR spectra were measured using an AVANCE600 spectrometer (Bruker BioSpin, Billerica, MA, USA) at a probe temperature of 10 • C. In this study, 1D imino proton spectra and 2D nuclear Overhauser effect spectroscopy (NOESY) spectra (mixing time of 150 ms were recorded using the jump-and-return or 3-9-19 pulse schemes for water suppression.

Native PAGE
The conformation of iMyo01-DNA, iMyo01-CS, 12nt-Ctrl, and 12nt-Ctrl-CS was analyzed by 20% native PAGE with or without 5 mM KCl.A concentrated KCl solution was added to the gel and running TBE buffer.The gels were stained with SYBR Gold (Thermo Fisher Scientific).

Multicanonical Molecular Dynamics (McMD) Simulation
The two molecules of iMyo01-DNA in explicit solvent were simulated by the trivialtrajectory parallelization (TTP)-McMD method [7,16].The simulation system (Supplementary Figure S1) contained two iMyo01-DNA molecules, 10,252 water molecules, 23 K + ions, and one Cl − ion to be neutralized electrically.The initial structure of one iMyo01-DNA molecule was built as a DNA helix model by NAB in AmberTools 22 [22].The other iMyo01-DNA was copied from the first and was separated by 10 Å.The force field of amber OL15 [23] was applied for iMyo01-DNA.The TIP3P model [24] was used for water.The periodic boundary condition was used with a truncated octahedron box.The long-range interaction was calculated by particle-meshed Ewald with a 10 Å cutoff.
The TTP-McMD [25] was conducted to sample the equilibrated conformations at 310 K. Before TTP-McMD, the volume of the system was equilibrated by MD simulation under the constant NPT condition of 1 bar.The energy range of the multicanonical ensemble covered from 270 K to 600 K.A total of 100 trajectories were used, and 14 iterative runs were performed to achieve a flat distribution along the energy range.The production run was conducted for 100 ns in each trajectory (total 10 µs).The snapshots were saved every 200 ps.A total of 50,000 snapshots were sampled, and the 988 conformations at 310 K were obtained as the 310 K canonical ensemble by the reweighting method.
The representative iMyo01-DNA structures were taken from the centroid of structural clustering using AmberTools22 [26], which was performed for the iMyo01-DNA molecules with the metrics of the root-mean square deviation of heavy atoms within iMyo01-DNA in the dimer state.A hydrogen bond analysis was performed using VMD tools version 1.9.4a53 [27].The structure images were generated by UCSF Chimera [28].

Statistical Analysis
The results were presented as mean ± standard error.Statistical comparisons were performed using multiple comparison tests with Dunnett's test following one-way analysis of variance.Statistical significance was set at p < 0.05.

Designing and Screening of iMyo-ODNs
The established 18-base myoDN, iSN04, has a tandem telomeric repeat (TTA GGG TGA GGG) as the core sequence for its myogenetic activity.In particular, the G-stack consisting of the 13-15th GGG is an indispensable structure for iSN04 activity [7].We designed eight non-telomeric PS-ODNs (iMyo-ODNs; iMyo01-iMyo08) with two G groups (GGG or GGGG) separated by a mini linker (T or TGA), thymines at the 5 ′ end (TT or TTT), and adenines at the 3 ′ end (AA or AAA) (Table 1).It has been speculated that iMyo-ODNs form a condensed conformation through interactions between 5 ′ -terminal thymines and 3 ′ -terminal adenines and have G-stacks within the structure, similar to that of iSN04.
To investigate the myogenetic activities of iMyo-ODNs, undifferentiated myoblasts were treated with PS-ODNs for 48 h and immunostained for MHC, a terminal differentiation marker of skeletal muscle.The signal intensities of MHC or the ratio of MHC + cells were automatically quantified in an unbiased manner using the screening system (Figure 1).iSN04 was used as a positive control, and iSN14 and iSN45 served as negative controls.In hMBs, iMyo01 and iMyo03 significantly increased MHC signal intensities to the same extent as iSN04.Other iMyo-ODNs and negative controls did not induce myogenesis (Figure 1A).In mMBs, iMyo03 and iMyo04 markedly increased the ratio of MHC + cells (Figure 1B).In chMBs, iMyo01-iMyo04 exhibited significant abilities to promote myogenic differentiation (Figure 1C).In all myoblasts, iMyo05-iMyo08 did not induce myogenesis, indicating that the latter GGGG, not GGG, is essential for iMyo-ODNs to function as myoDNs.
The effects of iMyo01 and iMyo03 on hMBs were further investigated by highresolution imaging and qPCR (Figure 2).Immunostaining revealed that both iSN04 and iMyo01 significantly induced the differentiation of hMBs into MHC + myocytes and multinucleated myotubes.iMyo03 also markedly accelerated myotube formation (Figure 2A).The mRNA levels of MyoD (MYOD1), a master regulator of myogenesis, were not altered by any of the PS-ODNs, but those of myogenin (MYOG), a myogenic transcription factor, were significantly increased by iMyo03 and were similarly induced by iSN04 and iMyo01.Transcriptions of myomaker (MYMK) and myomixer (MYMX), which are myogenin-regulated muscle-specific membrane proteins for myotube formation, were markedly upregulated by iSN04, iMyo01, and iMyo03 (Figure 2B).These results demonstrate that iMyo01 and iMyo03 as well as iSN04 promote the myogenic differentiation of hMBs.

iMyo01 and iMyo03 Bind to Nucleolin but Not to Berberine
Since iSN04 physically interacts with nucleolin to exert its myogenetic activity [7], iMyo01 and iMyo03 were tested for their ability to bind to nucleolin by precipitation assay.Biotin-conjugated PS-ODNs were immobilized on streptavidin beads.The soluble lysate of the C2C12 murine myoblast cell line was pre-pulled down with iSN14 beads to eliminate the absorption of non-specific proteins.After removing off-targets, the lysate was precipitated with iSN04, iMyo01, or iMyo03 beads.The precipitates were subjected to Western blotting to detect nucleolin.As shown in Figure 3A, nucleolin was precipitated by iMyo01 and iMyo03 as well as by iSN04.

iMyo01 and iMyo03 Bind to Nucleolin but Not to Berberine
Since iSN04 physically interacts with nucleolin to exert its myogenetic activity [7], iMyo01 and iMyo03 were tested for their ability to bind to nucleolin by precipitation assay.Biotin-conjugated PS-ODNs were immobilized on streptavidin beads.The soluble lysate of the C2C12 murine myoblast cell line was pre-pulled down with iSN14 beads to eliminate the absorption of non-specific proteins.After removing off-targets, the lysate was precipitated with iSN04, iMyo01, or iMyo03 beads.The precipitates were subjected to Western blotting to detect nucleolin.As shown in Figure 3A, nucleolin was precipitated by iMyo01 and iMyo03 as well as by iSN04.
iSN04 also interacts with berberine via its G-stack (GGG at the 13-15th base) in the presence of Ca 2+ or Mg 2+ [7].Berberine is an isoquinoline that binds to the G-quartet and stabilizes the G-quadruplex structure derived from telomeric sequences [29].The iSN04berberine complex shows higher myogenetic activity compared to iSN04 alone because berberine stabilizes the core G-stack within iSN04 [7].To determine whether iMyo-ODNs iSN04 also interacts with berberine via its G-stack (GGG at the 13-15th base) in the presence of Ca 2+ or Mg 2+ [7].Berberine is an isoquinoline that binds to the G-quartet and stabilizes the G-quadruplex structure derived from telomeric sequences [29].The iSN04berberine complex shows higher myogenetic activity compared to iSN04 alone because berberine stabilizes the core G-stack within iSN04 [7].To determine whether iMyo-ODNs form a complex with berberine, PS-ODNs and berberine were mixed in Ham's F10 and subjected to agarose gel electrophoresis.As shown in Figure 3B, the yellow emission from berberine (arrowhead) was detected at the same position as the red emission from iSN04 stained with EtBr (arrow), indicating the iSN04-berberine complex.However, all iMyo-ODNs, even the myogenetic iMyo01 and iMyo03, did not interact with berberine.This result suggests that the G-repeats of iMyo-ODNs form a different structure than the G-stack within iSN04.Since the myogenetic activities of iMyo01 and iMyo03 undoubtedly depend on their ability to bind to nucleolin, their conformations need to be determined directly.
stained with EtBr (arrow), indicating the iSN04-berberine complex.However, all iMyo-ODNs, even the myogenetic iMyo01 and iMyo03, did not interact with berberine.This result suggests that the G-repeats of iMyo-ODNs form a different structure than the Gstack within iSN04.Since the myogenetic activities of iMyo01 and iMyo03 undoubtedly depend on their ability to bind to nucleolin, their conformations need to be determined directly.

iMyo01 Forms a G-Quadruplex Structure
The imino proton NMR spectra of iMyo01-DNA were measured to characterize its structure (Figure 4A).In the absence of KCl, seven sharp (arrows) and four broad (arrowheads) signals were observed.After titration with 50 mM KCl, the NMR spectrum changed significantly, but seven sharp (arrows) and four broad (arrowheads) signals were observed, similar to those in the absence of KCl.The imino proton resonances between 10.5 and 12.5 ppm are characteristic for G-quadruplex formation [30], consistent with iMyo01 having no cytosine and therefore no G:C base pair.Since the G-quadruplex structure is stabilized by K + [31], it is considered that the KCl-induced shift of the imino proton spectrum of iMyo01-DNA represents the stabilization of its G-quadruplex by K + binding.To confirm G-quadruplex structure, the NOESY spectrum of iMyo01-DNA was measured (Figure 4B).The strong NOE signals between imino protons and NOE signals between H8 and imino protons strongly suggested that iMyo01-DNA forms the G-quadruplex (Figure 4B,C).

iMyo01 Forms a G-Quadruplex Structure
The imino proton NMR spectra of iMyo01-DNA were measured to characterize its structure (Figure 4A).In the absence of KCl, seven sharp (arrows) and four broad (arrowheads) signals were observed.After titration with 50 mM KCl, the NMR spectrum changed significantly, but seven sharp (arrows) and four broad (arrowheads) signals were observed, similar to those in the absence of KCl.The imino proton resonances between 10.5 and 12.5 ppm are characteristic for G-quadruplex formation [30], consistent with iMyo01 having no cytosine and therefore no G:C base pair.Since the G-quadruplex structure is stabilized by K + [31], it is considered that the KCl-induced shift of the imino proton spectrum of iMyo01-DNA represents the stabilization of its G-quadruplex by K + binding.To confirm G-quadruplex structure, the NOESY spectrum of iMyo01-DNA was measured (Figure 4B).The strong NOE signals between imino protons and NOE signals between H8 and imino protons strongly suggested that iMyo01-DNA forms the G-quadruplex (Figure 4B,C).

iMyo01 Forms a Homodimer in the Presence of K +
The NMR results showed the G-quadruplex structure of iMyo01-DNA; however, the 12-base iMyo01 sequence with only two G-repeats is too short to form a G-quadruplex, which requires at least 16 bases including four G-repeats [32].Therefore, it is assumed that iMyo01 forms a dimer to realize a G-quadruplex.The conformation of iMyo01-DNA was analyzed by native PAGE.As shown in Figure 5, iMyo01-CS (lane 6) was not stained with SYBR Gold due to its low number of purines (only three adenines), but other nucleotides were detected.iMyo01-DNA (lane 4) was observed in monomer mobility in the absence of KCl but was detected in dimer mobility in the presence of KCl, which is confirmed by the duplex of iMyo01-DNA and iMyo01-CS (lane 5).In contrast, no dimer mobility was observed for 12nt-Ctrl (lane 1), which never forms a G-quadruplex, even in the presence of KCl.The K + -dependent dimer formation strongly suggests that iMyo01-DNA forms a G-quadruplex structure in the presence of K + , in agreement with the results of the NMR.
Dimer formation of iMyo01-DNA was computationally simulated by McMD with the all-atom model.Using the simulated conformational ensemble, the distances between the centers of mass of the two iMyo01 molecules (d com ) were calculated.As shown in Supplementary Figure S2, their probability distributions showed that the dimer state (d com < 25 Å) is more stable than the separated monomer state (d com ≥ 25 Å).The conformation at d com ≈ 10 Å (arrow) had the highest probability, which was the most stable.Structural clustering of the dimer state represented the centroid structures of the top three clusters occupying 86.5% of all structures.Their d com values were 6.0 (42.1%, Figure 6A),

iMyo01 Forms a Homodimer in the Presence of K +
The NMR results showed the G-quadruplex structure of iMyo01-DNA; however, the 12-base iMyo01 sequence with only two G-repeats is too short to form a G-quadruplex, which requires at least 16 bases including four G-repeats [32].Therefore, it is assumed that iMyo01 forms a dimer to realize a G-quadruplex.The conformation of iMyo01-DNA was were detected.iMyo01-DNA (lane 4) was observed in monomer mobility in the absence of KCl but was detected in dimer mobility in the presence of KCl, which is confirmed by the duplex of iMyo01-DNA and iMyo01-CS (lane 5).In contrast, no dimer mobility was observed for 12nt-Ctrl (lane 1), which never forms a G-quadruplex, even in the presence of KCl.The K + -dependent dimer formation strongly suggests that iMyo01-DNA forms a G-quadruplex structure in the presence of K + , in agreement with the results of the NMR.6C).Since the conformation at dcom ≈ 10 Å showed the highest probability, those in the second cluster at dcom = 11.4 can be the stable dimer.
The analysis of the intra/intermolecular hydrogen bonds between the residues of iMyo01-DNA showed that the guanines form hydrogen bonds.The second cluster conformation indicated intermolecular hydrogen bonding between the G-repeats (Figure 6B), which may correspond to the G-quadruplex detected by NMR.On the other hand, the monomer state showed the intramolecular hydrogen bonds formed between G3-G5 and G7-G10 (Supplementary Figure S3), indicating that the monomer state had the stacking tandem guanines.Furthermore, the K + ions were found at the interface of iMyo01-DNA (Figure 6, purple spheres), which could reduce the electrostatic repulsion between them and thus contribute to the formation of the dimer.However, the G-stack in the dimer was slightly perturbed so that the simulated conformations did not have a fully coplanar arrangement of four guanines, suggesting that a more accurate force field of nucleic acids and solvents would be required to reproduce the complete G-tetrad.The analysis of the intra/intermolecular hydrogen bonds between the residues of iMyo01-DNA showed that the guanines form hydrogen bonds.The second cluster conformation indicated intermolecular hydrogen bonding between the G-repeats (Figure 6B), which may correspond to the G-quadruplex detected by NMR.On the other hand, the monomer state showed the intramolecular hydrogen bonds formed between G3-G5 and G7-G10 (Supplementary Figure S3), indicating that the monomer state had the stacking tandem guanines.Furthermore, the K + ions were found at the interface of iMyo01-DNA (Figure 6, purple spheres), which could reduce the electrostatic repulsion between them and thus contribute to the formation of the dimer.However, the G-stack in the dimer was slightly perturbed so that the simulated conformations did not have a fully coplanar arrangement of four guanines, suggesting that a more accurate force field of nucleic acids and solvents would be required to reproduce the complete G-tetrad.

Discussion
The present study successfully developed novel shorter myoDNs that promote myogenic differentiation.iMyo03 induced the differentiation of hMBs, mMBs, and chMBs as well as the established myoDN, iSN04, did.iMyo01 activated hMBs and chMBs, whereas iMyo02 and iMyo04 affected only chMBs.Such species-specificity of iMyo-ODNs may be

Discussion
The present study successfully developed novel shorter myoDNs that promote myogenic differentiation.iMyo03 induced the differentiation of hMBs, mMBs, and chMBs as well as the established myoDN, iSN04, did.iMyo01 activated hMBs and chMBs, whereas iMyo02 and iMyo04 affected only chMBs.Such species-specificity of iMyo-ODNs may be due to the differences in nucleolin structures.Although myoDNs are considered to interact with four RNA-binding domains (RBD1-RBD4) of nucleolin [9], the amino acid sequences of the nucleolin RBDs are not completely conserved among animals (Table 2).A previous study has also reported that G-quadruplex anti-thrombin DNA aptamers (15-TBA, 31-TBA, and RA-36) individually exhibited different activities against various mammalian thrombins [33].Clinically, aptamer cross-reactivity with the target ortholog in a patient can lead to side effects [34].Our results suggest that modification of the existing aptamer can alter not only the affinity but also the cross-reactivity to the target molecule, which needs to be considered in aptamer design.In the previous study, nucleolin inhibition by iSN04 increased p53 protein levels [7] because nucleolin binds to p53 mRNA and represses its translation [35,36].The present study confirmed that at least iMyo01 and iMyo03 physically interact with nucleolin.Thus, iMyo01 would increase p53 protein levels, resulting in the promotion of myogenesis because p53 cooperates with MyoD and the retinoblastoma protein to activate myogenic gene expression [37][38][39].Furthermore, nucleolin has been reported to target several mRNAs, including Bcl-2, interleukin 2, and matrix metalloproteinase 9, in addition to p53, in cancer [40].However, the role of nucleolin in myoblasts is poorly understood.Only one paper reported that miR-34b is involved in myoblast proliferation and differentiation via modulation of nucleolin protein levels [41].To clarify the mechanism of action of iSN04 and iMyo01, the functions and targets of nucleolin in myoblasts need to be elucidated in further studies.
The 12-base iMyo01 is the shortest myoDN found and is two-thirds the length of the established 18-base iSN04.Both iMyo01 and iSN04 act as nucleolin-binding aptamers, and their lengths are much shorter than the average of known aptamers (51 bases) [42].This is an advantage of our myoDNs as nucleic acid drugs.Because ODNs are chemically synthesized by sequential coupling of phosphoramidite monomers, longer ODNs have several issues related to purity, cost, and scalability for industrial mass production [43].As biomolecules, longer ODNs have a greater chance of being phagocytized and crossrecognized [15], compromising drug efficacy and safety.Therefore, shortening aptamers while maintaining their activities is an important and beneficial challenge in the development of aptamers as nucleic acid drugs.In addition, myoDNs must be taken up into cells because their target protein, nucleolin, is located in the nuclei.However, iSN04 is incorporated into the cytoplasm without a carrier [7]; such gymnosis generally occurs with ODN treatments [14].Since myoDNs need to cross plasma, endosomal, and nuclear membranes to reach nucleolin in the nucleus, the reduction in molecular weight will also contribute to their intracellular uptake.
On the other hand, the myogenetic activity of iMyo01 on hMBs was almost cellphysiologically equivalent to that of iSN04.The working doses of iSN04 and iMyo01 need to be carefully evaluated in further studies because iMyo01 acts as a homodimer.Although the nucleolin-binding ability of iMyo01 appeared to be lower than that of iSN04 in the precipitation assay, this is probably due to the inability of iMyo01 immobilized on the beads to efficiently form homodimers.
Both experimental and computational results indicated that iMyo01-DNA forms a Gquadruplex structure by dimerization in the presence of K + .Unmodified DNA molecules were used for these cell-free structural analyses due to the absence of nucleases.The nuclease-resistant PS-ODNs were used for cellular experiments, but they are mixtures of enantiomers of phosphorothioates and are not suitable for conformational analysis.Although the formation and stabilization of tertiary structures of nucleic acids are possibly affected by nucleotide modifications, a UV thermodynamic analysis has shown that phosphorothioates do not have a significant effect on the stability of the G-quadruplex structure [44].Therefore, the conformation of unmodified iMyo01-DNA confirmed by NMR and molecular simulation would essentially match that of PS-iMyo01.It is still difficult to selectively synthesize the designated enantiomer of PS-ODNs and is also hard to determine its structure experimentally.In silico calculations of configured enantiomers can be a powerful approach to overcome this problem and improve the productivity of aptamer design.The identification of iMyo01 in this study demonstrated that combined research technologies are useful for functional shortening of aptamers.It will hopefully contribute to the development of nucleic acid drugs.
For further applications of iMyo01, its efficacy and safety need to be evaluated in a variety of systems, including animal models.Myogenesis is a complex multi-step process from quiescent satellite cells to multinucleated giant myofibers [1].Previous and current studies have shown that iSN04 and iMyo01 promote myogenic differentiation and myotube formation of myoblasts [7,[9][10][11][12].However, a decline in the myogenic capacity of myoblasts is one, but not the only, cause of age-and disease-related muscle loss [2][3][4][5].Decreased numbers of satellite cells [45] and myofiber atrophy due to an imbalance between protein synthesis and degradation [46] are also major contributors to muscle wasting.Therefore, it is important to test myoDNs on satellite cells and myofibers to determine the therapeutic strategy.Except for our myoDNs, there are no aptamers that accelerate myogenesis, and there are no approved drugs for muscle wasting syndromes such as sarcopenia [47] and cachexia [48].Anti-myostatin antibodies are expected to treat muscle loss, but one of the antibodies, LY2495655, did not improve muscle mass and function at least in pancreatic cancer [49].For future success in reversing muscle loss, nucleic acid aptamers may be an indispensable drug modality.This study provides a candidate, iMyo01, and the wet and dry methodology for developing myogenic aptamers.

Conclusions
This study developed the 12-base short dimeric myoDN, iMyo01, which serves as an anti-nucleolin aptamer to promote myoblast differentiation as well as the well-known 18base iSN04 can.iMyo01 dimerizes in the presence of K + to form a G-quadruplex structure, resulting in nucleolin-binding and myogenetic capabilities.This is a successful instance of shortening the aptamer while maintaining its bioactivity, providing useful insights for the development of aptamers as nucleic acid drugs.

Figure 3 .
Figure 3. iMyo01 and iMyo03 bind to nucleolin but not to berberine.(A) Representative Western blotting image of nucleolin in the soluble whole-cell lysates of C2C12 cells precipitated by iSN04, iMyo01, or iMyo03.Beads control, streptavidin-coated beads without ODN.(B) Representative images of agarose gel electrophoresis of iSN04 and iMyo-ODNs mixed with berberine in Ham's F10.Ber, berberine.An arrow and an arrowhead indicate red and yellow emissions, respectively.

Figure 3 .
Figure 3. iMyo01 and iMyo03 bind to nucleolin but not to berberine.(A) Representative Western blotting image of nucleolin in the soluble whole-cell lysates of C2C12 cells precipitated by iSN04, iMyo01, or iMyo03.Beads control, streptavidin-coated beads without ODN.(B) Representative images of agarose gel electrophoresis of iSN04 and iMyo-ODNs mixed with berberine in Ham's F10.Ber, berberine.An arrow and an arrowhead indicate red and yellow emissions, respectively.

Figure 4 .
Figure 4. NMR analysis of iMyo01-DNA.(A) Comparison of 1D imino proton spectra of iMyo01-DNA in the absence of KCl and in the presence of 50 mM KCl. Arrows and arrowheads indicate seven sharp and four broad signals, respectively.(B) 2D NOESY spectrum (mixing time 150 ms).The NOE signals between imino protons (blue box) and the NOE signals between imino protons and H8 (red box).(C) Schematic representation of G-quartet and NOEs between imino protons (blue arrows) and between imino proton and H8 (red arrows).

Figure 4 .
Figure 4. NMR analysis of iMyo01-DNA.(A) Comparison of 1D imino proton spectra of iMyo01-DNA in the absence of KCl and in the presence of 50 mM KCl. Arrows and arrowheads indicate seven sharp and four broad signals, respectively.(B) 2D NOESY spectrum (mixing time 150 ms).The NOE signals between imino protons (blue box) and the NOE signals between imino protons and H8 (red box).(C) Schematic representation of G-quartet and NOEs between imino protons (blue arrows) and between imino proton and H8 (red arrows).

Figure 5 .
Figure 5. Native PAGE of iMyo01-DNA.Lane 1, 12nt-Ctrl; lane 2, duplex of 12nt-Ctrl and 12nt-Ctrl-CS; lane 3, 12nt-Ctrl-CS; lane 4, iMyo01-DNA; lane 5 duplex of iMyo01-DNA and iMyo01-CS; lane 6, iMyo01-CS, which was not stained with SYBR Gold due to its low purine content.Dimer formation of iMyo01-DNA was computationally simulated by McMD with the all-atom model.Using the simulated conformational ensemble, the distances between the centers of mass of the two iMyo01 molecules (dcom) were calculated.As shown in Supplementary Figure S2, their probability distributions showed that the dimer state (dcom < 25 Å) is more stable than the separated monomer state (dcom ≥ 25 Å).The conformation at dcom ≈ 10 Å (arrow) had the highest probability, which was the most stable.Structural clustering of the dimer state represented the centroid structures of the top three clusters occupying 86.5% of all structures.Their dcom values were 6.0 (42.1%, Figure 6A), 11.4 (27.0%, Figure 6B), and 18.6 Å (17.5%, Figure6C).Since the conformation at dcom ≈ 10 Å showed the highest probability, those in the second cluster at dcom = 11.4 can be the stable dimer.The analysis of the intra/intermolecular hydrogen bonds between the residues of iMyo01-DNA showed that the guanines form hydrogen bonds.The second cluster conformation indicated intermolecular hydrogen bonding between the G-repeats (Figure6B), which may correspond to the G-quadruplex detected by NMR.On the other hand, the monomer state showed the intramolecular hydrogen bonds formed between G3-G5 and G7-G10 (Supplementary FigureS3), indicating that the monomer state had the stacking tandem guanines.Furthermore, the K + ions were found at the interface of iMyo01-DNA (Figure6, purple spheres), which could reduce the electrostatic repulsion between them and thus contribute to the formation of the dimer.However, the G-stack in the dimer was slightly perturbed so that the simulated conformations did not have a fully coplanar arrangement of four guanines, suggesting that a more accurate force field of nucleic acids and solvents would be required to reproduce the complete G-tetrad.

BioTech 2024, 5 , 16 Figure 6 .
Figure 6.The representative iMyo01-DNA structures in the dimer state from McMD simulation.(A) The first cluster (42.1%; dcom, 6.0 Å). (B) The second cluster (27.0%; dcom, 11.4 Å). (C) The third cluster (17.5%; dcom, 18.6 Å).The clusters are shown as a rainbow-colored stick model from the 5′ to the 3′ end.The K + ions are shown as purple spheres.The hydrogen bond patterns in each cluster are plotted as a monochrome map.The gray scale indicates the probability of hydrogen bond formation between residues.The two iMyo01-DNA sequences are displayed sequentially on the axes.

Figure 6 .
Figure 6.The representative iMyo01-DNA structures in the dimer state from McMD simulation.(A) The first cluster (42.1%; d com , 6.0 Å). (B) The second cluster (27.0%; d com , 11.4 Å). (C) The third cluster (17.5%; d com , 18.6 Å).The clusters are shown as a rainbow-colored stick model from the 5 ′ to the 3 ′ end.The K + ions are shown as purple spheres.The hydrogen bond patterns in each cluster are plotted as a monochrome map.The gray scale indicates the probability of hydrogen bond formation between residues.The two iMyo01-DNA sequences are displayed sequentially on the axes.
. and T.T. are the inventors of Japanese Patent No. 7386507 covering iMyo-induced myogenic differentiation.Supplementary Materials:The following supporting information can be downloaded at https: //www.mdpi.com/article/10.3390/biotech13020011/s1:FigureS1:The simulation system of the iMyo01-DNA dimer; FigureS2: The probability of the center-of-mass distances between the two iMyo01-DNA molecules; FigureS3: The representative iMyo01-DNA structure in the separated monomer state.